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We expand previous work on an inverse approach to Einstein Field Equations where we include 
fluids with energy flux and consider the vanishing of the anisotropic stress tensor. We consider 
the approach using warped product spacetimes of class B\. Although restricted, these spacetimes 
include many exact solutions of interest to compact object studies and to cosmological models 
studies. The question explored here is as follows: given a spacetime metric, what fluid flow (timelike 

l""""" ■ congruence), if any, could generate the spacetime via Einstein's equations. We calculate the flow 

from the condition of a vanishing anisotropic stress tensor and give results in terms of the metric 
functions in the three canonical types of coordinates. A condition for perfect fluid sources is also 
provided. The framework developed is algorithmic and suited for the study and validation of exact 
' solutions using computer algebra systems. The framework can be applied to solutions in comoving 

and non-comoving frames of reference, and examples in different types of coordinates are worked 

Q ! out. 

O ' PACS numbers: 04.20.Cv, 04.20.Jb, 04.40.Dg 

m ; 

I. INTRODUCTION 

cr 

5^ ' The Einstein Field Equations (EFEs) have provided us with a wide variety of exact solutions that found application 
in the study and the understanding of black holes, compact objects models for neutron stars [jj], and cosmological 
models [H, Q . There are over a thousand of exact solutions to EFEs [l| including over 300 cosmological models [|J • 
, As discussed in [l[ and 0] , there is a need for studying the properties and interpretations of a large number of these 

J> ■ exact solutions. 

t-H \ In this paper, we expand the previous study by Ishak and Lake for fluids with zero energy flux. The framework 
developed there can be used to study current exact solutions, see e.g. [!, H|, or generate new ones, see e.g. Q. In this 
paper, we expand the work to fluids with energy flux and use the general form of the anisotropic stress tensor. 
Exact solutions with source fluids with isotropic pressure and non-zero flux have been used to model polytropic 
C"""- ' radiating stars and their gravitational collapse, see e.g. [3, El 0] • In cosmology, isotropic fluids with non-zero flux 
have been considered to address some problems in the the very early stages of the universe where dissipative processes 
were important, e.g. @, [l(| [HI- Such fluids were also considered, see for example [l2|, in order to model cosmic 
string networks with heat flux in inhomogeneous models. 

The usual approach to derive exact solutions is to start with a given spacetime metric, to specify a simple energy 
momentum tensor, often a perfect fluid, and to chose a comoving system of coordinates. These and other symmetry 
assumptions are usually necessary because of the difficulty in solving EFEs. The metric functions are then usually 
solved explicitly and then put back into the line element. We explore here an inverse procedure where some of the 
assumptions are relaxed. For example, we don't assume a perfect fluid source and we don't require a comoving 
coordinates system. For a given spacetime {M. , g) with manifold M. and lorentzian metric g, the inverse approach 
can be formulated as follows. What fluid flow, if any, could generate the spacetime (A^g) via Einstein's Equations? 

We explore this procedure using the warped product spacetimes of class B\ [ID, 0, [l5[. These spacetimes contain 
all spherical, planar, and hyperbolic spacetimes and all static spacetimes. For example, they include all Friedmann- 
Lemaitre-Robertson- Walker (FLRW) models, over a hundred cosmological models that are more complex than the 
FLRW models , neutron star relativistic models such as the Tolman Type- VII and the Buchdahl solutions [HI, [I?], 
H, El, H3| ! and large number of classical exact solutions such as the Schwarzschild and Reissner-Nordstrom black 
holes, the (anti) de Sitter spaces and many other spacetimes 

In the approach that we follow, the fluid velocity-field is derived from the metric subject to conditions on the fluid, 
rather than being put in by hand. In particular, the approach allows us to include solutions in either a comoving or 
non-comoving system of coordinates. Although the majority of exact solutions have been derived in a comoving frame 
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[l[ , solutions in the non-comoving frame have been been discussed in [2l|, [2^, [H, [24], [2^, [26| and were shown to often 
have rich kinematics, with shear, acceleration and expansion In this work, we are not particularly interested in a 
given type of coordinate system, however, we consider it an important advantage that the framework developed can 
handle exact solutions in comoving or non-comoving frames. 

Finally, this work provides an algorithmic framework to study exact solutions that is also suited for computer 
algebra system dedicated to exact solution studies like the online interactive general relativity database (GRDB) [27J • 



II. RELATIVISTIC FLUIDS AND WARPED PRODUCT SPACETIMES OF CLASS Si 

We consider in this paper, spacetime source fluids for which the energy momentum tensor can take the following 
decomposition relative to the fluid flow (velocity field) u a : 

Tab = (P + P)U a U b + pg a b + q a Ub + U a Qb + II a 6 (1) 

where 

p = T ab u a u b (2) 

is the energy density relative to u a , 

P=\{Tabh ah ) 1 (3) 

is the isotropic pressure relative to u a , and where h a b = u a u + g ao . 

q a = -T bc u b h ca , (4) 
is the momentum density which is also the energy flux relative to the fluid flow and satisfying q a u a — 0, and 

n a h = h^h^Td - -h ab h cd T cd (5) 

is the trace-free anisotropic pressure (stress)tensor satisfying H a a = 0, n a h = Il( a b), and H a bU b = 0. 

When q a — and n a h = 0, the energy momentum tensor T ab takes the usual form of a perfect fluid. 

The flow of the fluid source is specified by a congruence of unit time- like vectors u a = [u , u 2 ,0, 0] and is not 
necessarily expressed in the comoving system of coordinates. 

Next, we consider Warped Product spacetimes of Class B\ [H, 0, HE]. These spacetimes can be written in the 
form: 

ds 2 =dsl i (x 1 1 x 2 )+M 2 (x 1 ,x 2 )W 2 (x 3 ,x 4 )dsl 2 (x 3 ,x 4 ) (6) 

where signature(Si) = 0, and signature(S2) = ±2. Although very special, this class of spacetimes contains all of 
the spherical, planar, and hyperbolic spacetimes, including many spacetimes of interest to compact object studies and 
cosmological studies. 

We can further specify Si by writing: 

ds 2 ^ = A(x\x 2 )(dx 1 ) 2 + 2B(x 1 ,x 2 )dx 1 dx 2 + C(x\x 2 )(dx 2 ) 2 (7) 
where a,b, and c are functions of x 1 and x 2 only and the full metric is written as 

ds 2 = A(x 1 ,x 2 )(dx 1 ) 2 + 2B{x\x 2 )dx 1 dx 2 + C{x\ x 2 )(dx 2 ) 2 + M 2 (x 1 ,x 2 )W 2 {x3,x4:)((dx 3 ) 2 + (dx 4 ) 2 ) (8) 
The timelike condition on u a is 

u a u a =A{x 1 1 x 2 ){u 1 ) 2 + 2B(x 1 1 x 2 ){u 1 u 2 ) + C{x 1 1 x 2 )(u 2 ) 2 = -1 (9) 
and the spacetime (M, g) is time orientated by the restriction 

u 1 > (10) 
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Now, by the Einstein's Field Equations, we can write the anistropic stress tensor as 

n a f> = -h ( c h b] d G c d — -z-h ab h cd G c d (11) 

K iK, 

where G a b is the Einstein tensor, and its structure for the metric (JSj) is described in Appendix I, and where K = ^f- 
in full units. Note that for the spacetimes considered G33 = G44. 

In what follows, we will show that for the spacetimes under consideration, the velocity field can be determined from 
the vanishing of the anisotropic stress tensor (fTTj) . the timelike condition ([5]), and the restriction (JTDJ) . 

For the metric ([8]) and velocity field u a = [u 1 , u 2 , 0, 0], the nonzero components of U a b are Iln, II12, II22, n 33 , and 
II44 with II33 = II44, 

? 9 
II 1 H 

n n = (-) n 22 , n 12 = -(-)n 22 , (12) 



and 



where 



{Aiu 1 ) 2 + 2Bu l u 2 + C{u 2 ) 2 )M 2 W 2 
1133 - 2(ui) 2 (AC - B 2 ) 1122 (13) 



II22 = f(g ab ,Gb=g ac G cb ,u a ) 

2(u 1 f{AC - B 2 )((A(G 2 2 - G|) - BGDiu 1 ) 2 + {B(G\ + G| - 1G%) - AG\ - CG\)u x u 2 + {C(G\ - G|) - BG\){u 2 ) 2 ) 

3k(A(u 1 ) 2 - 2Bn 1 u 2 - G(« 2 ) 2 ) 2 

(14) 

Following previous work [HHH], we chose to use the mixed Einstein tensor G£ to express our results. 

When a comoving system of coordinates is allowed, i.e. u 2 — 0, the components of n o {, reduce to Iln = II 12 = 
and the non-zero components become: 

2{B 2 -AC){BG 2 -A{G 2 -GD) 
1122 = 3^ ■ (15) 

and 

w 2 r 2 (BG 2 -A{G 2 -G%)) , 
n 33 = ^ — \ iA 2 3JI = n 44 (16) 

For the spacetimes under consideration, the components of 11^ are proportional one to the other and solving for the 
components of u a that make the common factor equal to zero assures the vanishing of the whole anistropic stress 
tensor. 



III. FLUIDS WITH ENERGY FLUX AND VANISHING ANISOTROPIC PRESSURE 



It is of interest to discuss some physical situations where fluids with energy flux and isotropic pressure can be used 
to discribe such systems. A commonly used case in astrophysics is a polytropic radiating star [a. H^. For example 
neutron stars are well modeled by polytropes with index about in the range between n=0.5 and n=l. Also, for example, 
reference [28[ discusses that viscous, heat conducting fluids are less restricive than perfect fluids because they have 
more parameters and allow one to tune the geometric constraints with the physical constraints such as barotropic 
equations of state and other fluid thermodynamical constraints. Another situation is discussed in reference pj, where 
the authors explore spherically symmetric exact solutions with a barotropic equation of state, that are shear-free, and 
have heat flux. The solutions there were used in order to model the gravitational collapse of a radiating star which 
dissipates energy in the form of heat flux and where the heat flows from the central hotter region to the overlying 
cooler parts. 

In a cosmological context, it is well known that disspitative phenomena play an important role in early times of 
the universe. Refs 

discussed shear-free inhomogcncous models with heat flux and claim that the models can 
be used to address some the problems of the very early universe such as the horizon problem and entropy generation. 
Also, the authors of reference [i"2| derived exact solutions for isotropic cosmic string cosmology with heat flux in 
Bianchi type III space-time. As the authors motivate, it is widely accepted that, broken symmetries in the very early 
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universe generated topological stable defects such as domain walls, strings, and monopoles [29[. Cosmic strings have 
attracted most interest among these and are still subject to searches with most recent observations (30l. l3l|. 

A physically reasonable energy-momentum tensor has to satisfy the dominant energy condition (DEC) [32]. This 
means that any observer with a timelike 4-velocity measures a non-negative local mass-energy density, the pressure 
must not exceed the mass-energy density, and the energy flow vector Q b = T ab u a be non-spacelike. For our models, 
these conditions arc satisfied if 

P > (17) 
P + P > (18) 
p-p > 0. (19) 
Q b Q b < (20) 

where p and p are given by @, ([3]), respectively. A less restrictive condition is the weak energy condition (WEC) 
requiring only inequalities (|17[) and (|18[) . 

We can see from (JU)-© that once we derive the fluid velocity field that gives a vanishing anisotropic stress tensor 
for the source of a spacetime with a known metric, we can calculate the remaining fluid properties. This will also 
allow us to express the strong and weak energy conditions in terms of functions of the metric: 

{G ab u a u b ) = Gn w V + G 22 u 2 u 2 + 2Gi2 uV > 0, (21) 
(G ab u a u b +\G ab h ab ) = Gn u 1 u 1 +G 22 u 2 u 2 +2G 12 u 1 u 2 +l-(G 11 h n +G 22 h 22 +2G 12 h 12 +G 33 h 33 +G M h 44 ) > 0, (22) 

(G ab u a u b ~\G ab h ab ) = Gn wV+G 22 uV+2Gi 2 uV-i(Gn h 14 +G 22 h 22 +2G 12 h 12 +G 33 h 33 +G iA h 44 ) > (23) 
and 

nb n -2{G l2 v} + G 22 u 2 )(Gii U 1 + G 12 u 2 )B + A(G 12 u 4 + G 22 u 2 ) 2 + G(Gn U 1 + G 12 u 2 ) 2 

Q Qb ^ [AC — B 2 ) " ° (24) 

where Einstein's Field Equations have been used to give the results in terms of the Einstein tensor which can be 
calculated directly from the metric. An expanded discussion of the energy conditions in Warped Product type B\ 
spacetimes can be found in [l3|, 0, HB| • 

IV. DOUBLE-NULL COORDINATE FORMALISM 

We solve for the fluid flow and its phenomenology in the three canonical types of coordinates. We start in this section 
with the double-null coordinate representation (Kruskal-Szekeres type-like) characterized by A(x 1 , x 2 ) = C(x 4 , x 2 ) = 
into (J?]). This gives the metric the form: 

ds 2 = 2B{x\x 2 )(dx 1 dx 2 ) + M 2 {x 1 ,x 2 )W 2 (x3,xA)((dx 3 ) 2 + (dx 4 ) 2 ) (25) 

We can always choose bix 1 , x 2 ) < 0, and with the restriction (jTUJ) , it follows from the timelike condition that 

" 1 = - 2B ( ^y >ft (26) 

Note that the u 2 = 0, i.e. the comoving case, is not possible here. 

As mentioned earlier, each component of H ab is proportional to the others and as u 2 ^ we can choose the simplest 
component IIn, and set that equal to zero. We use ([26]) to substitute for u 1 and also note that, for this canonical 
type of coordinates, the Einstein tensor has G\ = G 2 , also noted in [3[. 

We find that the condition for iln = 0, and therefore also for H ab = 0, is as follows (again, following previous work 
[3, HI], we use the mixed Einstein tensor G b = g ac G cb to express our results): 

B 2 G\{u 2 ) 4 + (u 2 ) 2 (G{ - G%)B + jGl = (27) 
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We solve for u 2 and u 1 and, by (|26|) and (|10p . we keep only the positive solutions: 

(1/2) 



which gives us: 



2 / (gj ~ G\) + y^glg - 2G\Gl + (G|) 2 _ G i G 2) 
" - 2^ (28) 



/ G 1 \ (1/2) 

u i = U2 (29) 

I 2B[{G\ - G\) + ^((G}) 2 - 1G\Gl + (G3) 2 - G*G 2 )] i 



(1/2) 



2 _ I (gj - gj) ~ V((Gj) 2 - gGjgl + (G3) 2 - G\G'{) . 



which gives us: 



/ G 1 \ (1/2) 

u 1 = ° 2 (31) 

V2S[(Gi - G\) - V((G}) 2 - 2G1G| + (G|) 2 - G^G 2 )] ^ 

For the metric (|25|) , the fluid flow u a given by equations ([28l [29]) and ([30l l3T|) assures the vanishing of the anisotropic 
stress tensor. 

The next question of interest is what further geometrical condition, in terms of the metric functions, is required in 
order for the fluid to be perfect. Thus, we are interested in finding the sufficient condition to make the fluid flux-free. 
By the EFE and by (|3"2"|) , we can write the condition in terms of G^ and u a as 

q a = --h c a G cd u d = 0. (32) 

K 

For the spacetimes considered q a has the form q a = [q 1} q 2 , 0, 0]. From (|3"2"|) and after removing the square roots by 
multiplying each component by its radical conjugate, we can solve for the geometric condition that makes the flux 
equal to zero. Comparing the q\ and q 2 components, we find that the condition for zero flux is: 

G\(-G\G\ + {G\f - G\G\ - G\Gl + G\G%) = (33) 

which we can rewrite as 

Gl(EECR) = (34) 
where, EECR is a closure relation for the Einstein Field equations, and reads 



EECR = G\G\-{G\-Gl){Gl-Gl) (35) 

So, with u a given by (|2"51 12"9"| or ((501 |3"T) (thus assuring a vanishing anisotropic stress tensor), it turns out that 
the condition EECR — is the necessary and sufficient condition for a zero energy flux fluid and thus a perfect 
fluid. This result complements previous work by [33| for spherical symmetric spacetimes and that of reference Q for 
warped product spacetimes. Interestingly, in the approach followed in [3j, the velocity field was first solved using the 
zero energy flux condition and then the same perfect fluid closure condition was found to assure the isotropy of the 
pressure. It becomes clear then that the EECR is a closure condition for the EFEs solved with a perfect fluid energy 
momentum tensor for the spacetimes under consideration. 

As mentioned earlier, one of the goals of developing such an inverse-approach framework is to study and validate 
existing exact solutions. Let's apply it here to the following exact solution in double-null coordinate systems. 

The Einstein-de Sitter universe in double null coordinates fs^ j 

ds 2 M = C 2 (u + v) 4 (-dudv + ^~^ 2 rf» 2 ), (36) 
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where C is a constant and <Kl 2 is the metric of a unit sphere. 

In this spacetime B{x l ,x 2 ) = - c ' 2 ^+ v )\ M(x 1 ,x 2 ) = c("+^ 2 ("-^) ; W(x3,xA) = 1, C is a constant and dft 2 is the 
metric of a unit sphere. Now, Using (j28[) and (129)) and the metric, we derive algorithmically the flow that assures 
n a (, = and find u 2 = c ^ +v yi and u 1 = c ^ +v y2 ■ This flow coincides exactly with the solution derived in the 

original reference [34J. Further, we verified that the closure relation is satisfied, i.e. EECR — 0, so the source is 
flux-free and thus is a perfect fluid as it should. 



V. NULL COORDINATE FORMALISM 

We consider here the analysis using the canonical null coordinate system (Bondi-like coordinate systems) with 
C(x 1 ,c2) = or A(x 1 ,x 2 ) = 0. We choose C{x 1 ,x 2 ) — and A{x 1 ,x 2 ) < 0), so the metric for a Warped Product 
Type B\ spacetime takes the form: 

ds 2 = A(x 1 ,a^)(dx 1 ) 2 + 2B(x\x 2 )dx 1 dx 2 + M 2 (x\x 2 )W 2 (x3,xA)((dx 3 ) 2 + (dx 4 ) 2 ) (37) 

As also noted in [3j, the components of the Einstein tensor are related here by 

B(x 1 ,x 2 )(G\ — Gf) = A{x 1 : x 2 )G\. (38) 

In order to solve for u 1 and u 2 , we substitute the timelike condition ([9]) and equation (|38[) into n o {, = 0. This time, 
however, we use the component II22 as it is not directly proportional to u 2 . The condition for II22 = 0, and thus 
Ii a b = is found to be: 

4B 4 G?(« 2 ) 4 G^ + 2{{u 2 ) 2 G\{G\ - G%)A + Gi{G\ - 2G% + G^))(u 2 ) 2 B 3 + {-{ u 2 f{G\f A 2 - 4(G^G? - (l/4)(Gj - G 2 2 ) 
(G\ - 2G\ + G 2 2 ))(u 2 ) 2 A + (G?) 2 )B 2 + 2(~(u 2 ) 2 G 1 2 (G{ - G%)A + G?(G| - G%))AB + A 2 (A(u 2 f(G 1 2 ) 2 + (G| - G 2 2 ) 2 ) = 0(39) 

Solving (f39|) and using the timelike condition ([9]), we find the solutions ([40|) - (|46]) below for the flow. For each u 2 
found, solving the timelike condition © for u 1 gives two possible solutions, however, we note that only one of them 
satisfies the restriction (fTOj) . For brevity, We write here only the two positive solutions (i.e. with u 2 > 0). 

X = 4G 2 (G| - Gj)B 3 + 2b 2 G\G\A + 2G\A 2 (Gl - G\)B (40) 
Y = 16(G 2 ) 2 B 5 (((G 2 -G\) 2 - G 2 G\)B + (G 2 - (%)AG\) (41) 
Z = 2B 2 G 1 2 (4B 2 G 2 1 + A 2 G\) (42) 



^ - ,M±£ (43) 



and 



and the second set: 



and 



2 x-Vy 

u =y^z— (45) 



(46) 
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We verified that all our solutions for u a do give back H a b = 0. We also found from (fTSf . (fTBf . (|43f and (|45f that a 
comoving coordinates system w 2 = is possible if and only if 

BG\ - A(G\ - G%) = 0. (47) 

Now, let us go back to the more general solutions (|4"3")) - (|4"5")) above for the flow, and consider further conditions for 
the fluid to be perfect. Again, with vanishing anisotropic stress tensor, we are interested in the zero-flux condition. 
We write out q a and substitute our solutions for u 1 and u 2 into the results and multiply by the radical conjugates to 
get rid of the radicals so we can search for the geometric condition to make q a — and we find the condition to be 

{EECR){{{G\fGlGl + {G\) 2 G\Gl + {G\fG\ - {G\f{Glf - {G\) 2 G\Gl)A 6 + {AG 2 Gl(G\) 2 - A(Gl) 2 G 2 G 2 )BA 5 + 
( - AG 2 G\{Gl) 2 + 8(G 2 ) 2 (Gl) 2 + %G 2 G\G 2 Gl - AG 2 G\(G 2 ) 2 ^B 2 A i + \Q{G 2 fB i A 2 G\ + 16(G 2 ) 4 B 6 ) = 0, (48) 

where EECR is as given in equation (|35|) . So, again we see that the closure condition, EECR, if zero, will make the 
flux go to zero, but it is not the only factor that can be responsible for zero flux. Therefore, we can only show here 
that EECR = is a sufficient condition for anisotropic stress tensor. 

Now, let's apply these results to exact solutions and metrics in Null coordinate systems: 

The Davidson solution in null coordinates jizdi] 

ds 2 M = -2H(u, r)(du) 2 - 2dudr + ur 2n {(dx) 2 + (dy) 2 ), (49) 

where H(u,r) = + ^kr" 1 ^™- 1 ^ and n = m(m — l)/2. i.e., A(u, r) = —2H(u,r), B = —1, r(u,r) — y/ur n , and 
w(x, y) = 1. 

Using our equation (|40p - ([4^]) . we verified that the velocity field for which n a {, all reduce to u 2 — and u 1 = 1 
as claimed by the authors. Moreover, we verified that our equation (|4T|) is satisfied and therefore the anisotropic 
stress tensor must vanish in the comoving frame and is consistent with the explicit subsitutions. Furthermore, 
we verified that the EECR condition is satisfied so the source is flux free and thus a perfect fluid, in agreement with [2(| . 

The General Bondi metric |H/ 



ds 2 M = c 2 (v, r)f{v, r){dvf ± 2c(v, r)dvdr + r 2 dtt 2 (50) 

For this spacetime, A(t,r) = C(v,r) 2 f(v,r), B{t,r) — c(v,r), M(t,r) — r 2 , and W(9, 4>) — 1 and (+) is for 
advanced and (-) is for retarded null coordinate v. We find from the zero anisotropic stress tensor condition: 



u 2 = 



-(c, r -c)(3r 4 /, r c, r +2r A fc,rr +r 4 cf, r 



-Afr 3 c r 



-4/r 2 c - 



2c) 



8r 2 (c, r r — c) 2 



(1/2) 



(51) 



and 



! _ f -(c, r -c) (3r 4 /, r c, r +2r 4 /c, rr +r 4 c/, rr -4/r 3 c, r -4/ r 2 c - 2c) \ (1/2) 
U ' I 8r 2 (-c/) 2 (c, r r- C ) 2 ) ^ 

/ -(c, r -c)(3r 4 /, r c, r +2r 4 /c, rr +r 4 c/, rr -4/r 3 c, r -4/r 2 c - 2c) ' N (1/2) 



8r 2 (-cf) 2 (c, r r-c) 2 ' {-cff 

where a "," followed by a subscript coordinate denotes a partial derivative with respect to that coordinate. We 
verified that Il a b = for all our solutions to u a . However, when we substitute the comoving solution (u 2 — 0) into 
H a b we see that we get zero anisotropic pressure only if the metric satisfies the condition A(G\ — G 3 ) — BG\ = 0. 
The limitations of the comoving frame for perfect fluid sources of this metric have already been noted in [f| . 
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VI. DIAGONAL COORDINATES 
A. formalism 

The third canonical coordinate system that needs to be covered is the diagonal system where B(x ,x ) = and 
the metric takes the following form. 

ds 2 = A(x\x 2 )(dx 1 ) 2 + C{x\x 2 )(dx 2 ) 2 + M 2 (x\x 2 )W 2 (x3,x4)((dx 3 ) 2 + (da; 4 ) 2 ) (53) 

We can always choose C(x 1 , x 2 ) > and A(x x , x 2 ) < 0. As noted in the Einstein tensor components in diagonal 
coordinates satisfies 



We use the timelike condition and the relation above into II22 = in order to solve for u 1 and u 2 . The condition for 
II22 = 0, and thus for n a fc = 0, is: 



(1 + C{u 2 f)(C{u 2 ) 2 {G\ - G 2 ) - (G 2 - G%) + 2u 2 GW-A(l + C(u 2 ) 2 )) = (55) 

A lengthy but straightforward calculation gives the following pairs of solutions (we show only the positive solutions 
for brevity, i.e. u 2 > 0) 

(1/2) 



2 _ ( (gj - gjXgj - G\)C 2(G\) 2 A + VgGjgqgijg + C{GIG 2 + G\G\ - {G\) 2 - G 2 G\)) . 
U { {A{G\) 2 AC +{G\-G 2 ) 2 {C) 2 ) j ' 1 1 



1 _ ( 4( V(G 2 )M(-(G| - G 2 )(Gl - G\)C + (Gl) 2 A) {G\ Gl){G% G\)C + 2(G\) 2 A) \ U/2) 



u = 



(~A){2{G\-Gl) 2 C + %{G\) 2 A) 



(57) 



and 

2 _ ( {G\-Gl){Gl-G%)C-2{ G\fA - 2^ A{G\) 2 {{G\) 2 A + C{GlG 2 + G$G\ - (Gi) 2 - G 2 G{))\ ^ 



(4(G 2 )MC + (G{- G 2 ) 2 {C) 2 ) 



(58) 



(1/2) 



1 _ [ 4(- V(GWA-(G 3 3 - mm G\)C + {G\) 2 A) (G\ - G 2 ){Gl - G\)C + 2{G\) 2 A) ■ 

^ (-A)(2(G\~G 2 ) 2 C + 8(Gl) 2 A) ) [ ) 

Again, we verified that all the components of Tl a b vanish as they should when we substitute the solutions found. 

Just for completeness, let's verify the requirement for comoving coordinates using our more general formalism. This 
follows from solving for u 2 — in our general solutions ([5^)1 or l[58p and is found to be G\ = G\. 

Conversely, if one chooses the comoving frame with starting point u a = ( , 0, 0, 0) then II11 = II12 = II21 = 
and we are left with 

II22 = 2C(Gi " Gl) (60) 



and 



n ss = n u = - M2w 'V- Gi) 



as the nonzero components. We see that the necessary and sufficient condition for the vanishing anisotropic stress 
tensor is G 2 = G3, in accord with above. This also verify that the general framework developed here covers both 
non-comoving and comoving frames of references. 
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Now, the next step of interest is to solve for geometrical conditions that will make the flux go to zero when the 
anisotropic stress tensor is zero, thus giving a perfect fluid. As mentioned earlier, the flux vector q a has two components 
and we find that the condition for the vanishing of q\ and (72 is 

(£i5Ci?)((Gi-Gi) 2 (G^ 

(62) 

with EECR given by equation (155)) . 

We see once again that the closure condition, EECR, is a common factor, but not the only common factor between 
the components of q a . Thus we can only show here that EECR = is a sufficient condition for the fluid to be perfect 
when the velocity vector field is chosen to make the anisotropic pressure tensor zero. 



B. Application of the formalism to known exact solutions 

We apply our formalism to exact solutions in Diagonal coordinates. The formalism can handle solutions with or 
without energy flux, in comoving or non-comoving frames of reference. 

Sussman shear-free solutions with heat flux 

l2 —N 2 (t, r)dt 2 + dr 2 + r 2 (d0 2 + sm 2 (6)d(j) 2 ) , . 

ds = lWJ) (63) 

where 

N(t,r) = k^VM^WH-Wr 2 ) + fc2e -^(7(t)+|«(t)r 2 ) ) (64) 

and 

L(t,r) = ^eV^hW+l^^ 2 ) + fc 4e - y/^hW+^Wr') . (65 ) 

So A(t,r) = -j^yr , C(t,r) = -jjjjyj, M(t,r) = x^y, and W(9,<f>) = 1. Using our general equations |(5BJ)-||5SJ) for 
zero anisotropic stress tensor, we find that all the solutions for the components of the velocity field reduce to those of 
the comoving asynchronous frame u a = ( jm^k ; 0, 0, 0) and the solution has radial flux q a = [0, q r ,0, 0]. The EECR 
condition (|35| is not satisfied and thus is consistent with the non-zero flux found. All results in agreement with [281 ] . 

The Dadhich and Patel shear-free solution with heat flux fUl l 



ds 2 = -(r 2 + P{t)) 2n dt 2 + (r 2 + P(t)) 2m (dr 2 + r 2 (d8 2 + sm 2 (0)# 2 )) (66) 

This spacetime correspond to A(t,r) = -(r 2 + P(t)) 2n , C(t,r) = (r 2 + P(t)) 2m , M(t,r) = r(r 2 + P(t)) m , and 
W(9,(f>) — 1. We only considered their particular solution with 2m = 1— sqrt(3/2) < and 2n = sqrt(3/2) in the 
comoving frame. Using our general equations, we verified that the fluid has zero anisotropic stress tensor and radial 
flux. However, the general solution given there was found to have some problems (e.g. not shear- free and G\ ^ G3). 
We suspect a misprint in the second part of their equation (1), i.e. 2n — 2m ± sqrt(8m 2 + 8m + 1) as it does not 
reconcile with the working solution given by equation (2) there. 

The Govender and Govinder shear-free solution with heat flux fldt l 



ds 2 = -A 2 r Sa /( 1+a Ut 2 + (c t + Cl ) 4 / 3 ( 1+ «) 

Again using our formalism, we find that the general solutions (p)6"|) - (l5i?|) for zero anisotropic stress tensor give u 2 = 0, 
i.e. a comoving frame of reference, consistent with G\ = G3. The fluid has radial flux, consistent with a non-zero 
EECR condition. All in agreement with [UJ. 



4(1 



~dr 2 +r 2 dQ 2 



(67) 
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The FLRW metric in diagonal coordinates. 



ds 2 M = -{dtf + a ® { f r ] + a(tfr 2 dn 2 (68) 
1 — kr A 

2 ( t) 

In this spacetime A(t,r) = — 1, C(t,r) = i_ kl ! 2 , M(t,r) — a(t)r, and W(6,(p) — 1. It is a trivial exercise to verify 
that the source-fluid is a perfect fluid in the comoving frame. Our four solutions for u 2 in this spacetime reduce to 
the same solution: u 2 — 0, and the corresponding u 1 — 1, as expected. 

The McVittie- Wiltshire metric with e = +1 in non-comoving diagonal coordinates. Fiji / 



ds 2 M = e 2 ^ z \-{dn) 2 + (dO 2 + fdil 2 ) (69) 

where z = e(£ 2 — n 2 )/^, and /3(z) is an undetermined function of z with (3 ZZ — fi 2 ^ 0. Using the formalism developed, 
we find that the anisotropic pressure tensor vanishes for the following velocity field given in the non-comoving frame 
of reference by 

u 2 = — — — -p and u 1 — — — — ^ (70) 

which agree with the forms found in Also, we verified that the condition EECR is satisfied so the fluid is flux free 
and thus perfect in agreement the original reference (25j . 

The Senovilla-Vera G2-cosmological models in diagonal non-comoving coordinates \2_ 

ds 2 M = -{dtf + (d X ) 2 + C0S IZ^ X) J ^ 2 + C °t^<\^ 2 ^ 

cosh [fit) cosh (/it) 

We use our general solution (|56]) - ([59^) to determine the fluid flow for which the anisotropic stress tensor is zero and 
find 

u 2 = ( (1-C0s(^) 2 )c0sh(^) 2 \ (1/2) 

\ (2 cos((j,x) 2 cosh(fit) 2 — cos(fix) 2 — cosh(pt) 2 ) J 

and 

u i = ( (cosh( M t) 2 -l)cos( M x) 2 \ (1/2) 

\ (2 cos(fj,x) 2 cosh(/ii) 2 — cos(px) 2 — cosh(/it) 2 ) / 

which after transformations match exactly the solutions provided by the authors of [2l|, i.e. 

9 tan(ux) , , tanh(ut) .„,. 

u = = ' and u 1 = = ; - (74) 

•wtanh 2 (/it) — tan 2 (/iic) A/tanh 2 (/it) — tan 2 (/ix) 

Furthermore, we verified in this non-comoving frame that the condition EECR is satisfied so the fluid is flux-free and, 
thus, perfect in agreement with [2lJ. 



C. Recovering some solutions from application of the formalism 

As mentioned in the introduction, our formalism has been developed so it can be used to validate and study 
existing solutions, however, it could also be used, if desired, in order to derive new solutions. Such attempts have 
been explored for example in the previous papers [H, [B| . We give here an idea on how our formalism can be used to 
recover some solutions without the need to solve explicitly the Einstein Field equations. 
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Example-1: A perturbed Static Spherically Symmetric Line Element 

Heintzmann [3^ | derived a static spherically symmetric exact solution with a perfect fluid source that was shown 
in (37j to fit several physical acceptability conditions. We consider this solution and add a time perturbation to the 
metric via the functions R(t) and S(t) as they appear in the following line element 



( - — + + **** (a 2 (1 + ar 2 ) 3 + S(t))dt 2 . (75) 

r o 1 1 2 



ds 2 = 

\ 1 

2 1+ar 2 



The original static solution has isotropic pressure 



(7 ar 2 C + 3 Vl + 4 ar 2 ar 2 + C - 3 Vl + 4ar 2 ) 



P= ~ 3 / 2 ^ n 2 ,2 n A j — > (76) 

K (1 + ar 2 ) Vl + 4ar 2 



energy-density 



a ( (1 + 4 ar 2 ) 3/2 ar 2 + 9 ar 2 C + 3 (l + 4 ar 2 ) 3/2 + 3 c) 

P = 3/2 =7= = ^ , (77) 

(l + 4ar 2 ) 3/2 (l + ar 2 ) 2 K 

and zero energy flux. However, with the time perturbations above, the solutions becomes significantly complex with 
energy density given by 

p per = (24 R (tf a 3 r 6 C - 9 C 2 \J\ + 4 ar 2 a 2 r 4 R (t) + 24 R (t) + 4 ar 2 ar 2 + 4 R (t) \/l + 4 ar 2 

-18 a 2 r 4 R (t) C - 12 R (tf ar 2 C ~ 42 R (tf a 2 r 4 C + 12 R (tf ay/l + 4ar 2 r 2 
+4 R (tf \/l + 4ar 2 + 27 r 4 a 2 R (t) \/l + 4ar 2 + 18 ar 2 C ~ 20 r 6 a 3 R (t) \/l+4ar 2 
+4 R (tf a 3 r 6 y / l + 4ar 2 - 15 R (tf a 2 r 4 a/I + 4ar 2 + 54 a 2 r 4 C + 24 a 3 r 6 y/l + 4 (a) r 2 
+18 Vl + 4ar 2 ar 2 + 78 a 2 r i \ / l + 4ar 2 - 72 a 3 r 6 CR (t) + 9 V'l + 4ar 2 .R (t) 2 a 2 r 4 C 2 ) / 

( (-2 R(t)y / l + 4ar 2 + i? (i) \/l + 4ar 2 ar 2 + 3 R (t) ar 2 C - 2 \f\ + 4ar 2 - 2 y / l + 4ar 2 ar 2 ^j 2 y / l + 4ar 2 r 2 h 

(71 

and an anisotropic pressure with pi ^ p2 and both too cumbersome to display here in the paper. Moreover, the 
energy flux is now not zero and has one radial component given by 



q r = - (-2 VT+TaT 2 + VTTTaT^ar 2 + 3 ar 2 c) — R (t) j (k (-2 R (t) VT+TaT 2 + R (t) vT+TaT^ar 2 

+3 R (t) ar 2 C - 2 a/I + 4 ar 2 - 2 a/1 + 4 ar 2 ar 2 ) 2 r V^ 2 + 3 A 2 (a) r 2 + 3 A 2 a 2 r 4 + A 2 a 3 r 6 + S (tfj . (79) 

Now, using our formalism, we are interested in finding what conditions are required from the perturbed metric in order 
to have a perfect fluid source solution. In this case, i.e. diagonal and comoving frame, our formalism indicates that 
the condition ([35]) is sufficient to assure a zero flux source and gives the condition dR(t)/dt = 0, i.e. R = Const = Rq. 

Next, from our equation ([SO]) , the condition for isotropic pressure in this case is G\ — G\ = 0. Using R(t) = Ro and 
after some algebra, we factor this equation into a cumbersome but informative form as 

F(r) S(tf + G(r) S(t) + H(r) R = 0. (80) 

Since S(t) can be a function of t only, the only solution is given by S(t) — Const = So- This forces the solution back 
to a static one. The Heintzmann solution [3|| is a special case with So = Ro = 0. 

Example-2: Einstein-de-Sitter in non-comoving double-null coordinates 

We alter the metric derived in (34| and use the anzats r(x 1 , x 2 ) = C(Av + u)(u + v) 2 /2 giving the line element 

ds 2 M =C 2 (u + v) 4 (-dudv + ± U)2 aT> 2 ), (81) 
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where C is a constant and <Kl 2 is the metric of a unit sphere. 

With the flow vector field u 2 = c ^ +v yi and u 1 — c ^ +v yi , the fluid is not pressure or flux free and has non- vanishing 
anisotropic pressure. 

Using our formalism we will find what conditions are required in order to recover the Einstein-de-Sitter dust universe. 
First, without solving any Einstein differential equations, and directly from our derived equation (127j) for vanishing 
anisotropic pressure, we obtain the first constraint 

(A + 1) (2 Av + v + 3 u) (14 A 2 v 2 + 3 v 2 A + v 2 + 2 vuA 2 + 30 Auv + 4 uv + 15 u 2 + 3 u 2 A) 
-16- i-i n ^ = 82 

(Au + w) 4 (u + i;) u C 4 



and from our equation (|35p for zero flux, we obtain 

2 (^+l)(2^ + , + 3 U ) = 
(Aw + u) 2 (u + w)° C 2 

A closer look to the two equations above, and for arbitary coordinates u and v, it follows that both constraint 
are verified for A — —1 and a further check gives p — ■^{G a bh ah ) = and q a — —^Gb c u b h ca . This recovers the 
exact solution presented in (34|. The energy conditions reduce to p = ^G a bU a u b = 48/ n(u+v) 2 C 6 > 0, clearly satisfied. 

Example-3: Friedmann-Lemaitre-Robertson- Walter (FLRW) in comoving diagonal coordinates 

Here we alter the FLRW metric by replacing the function a(t) by a function f(t) in the metric components geg, 
and giving the line element 

ds 2 M = -{dtf + a( * r f/ + f(t) 2 r 2 dn 2 . (84) 

Now, in the comoving frame, u a — (1,0,0,0), we would like to recover the FLRW metric without solving the 
Einstein equations but rather by using some results of our formalism. The condition (|35p is sufficient to assure a zero 
flux source and gives the constraint 

-^/(*) + ^*(t) = 0. (85) 

The solution to ([55]) follows immediatly as 

/(*) = Const x a(t) (86) 

but is not enough to provide isotropic pressure. For that, we use the condition G\ = G3 from our equation (f6"0"l) which 
is sufficient to assure the vanishing of the anisotropic pressure and gives the constraint 

if H)) 2 -fit) p it) ~ r" 2 ) a it) + f (t) (|/ (t)) |a (t) + + £a (t)) (/ (t)) 2 = 0. (87) 

Putting equation (|86|) into the second constraint ([87]) yields 

a(i) 2 (Co?isf - l) 2 = 0, (88) 

giving thus the desired result, i.e. Const = 1 or f(t) = a(t) by (f86|) . This restores the well know FLRW metric with 
a perfect fluid source and energy conditions 

P = 3 L2 >0 (89) 

K (a it)) 

a(t)-Ca(t) + 2 (±a(t)) 2 + 2k , x 

p + p= 2 W ^ W Vdt 2 V y; > (90) 

K (a(t)) 

p-p = 2 Kdt K>> y g W > 0. (91) 

K (o(t)) 
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VII. SUMMARY 

This work complements previous results in [3J] where the flow was calculated from the zero energy flux condition 
and then perfect fluid conditions were derived. The current work starts with the derivation of the flow for which 
the anisotropic stress tensor vanishes while allowing for energy flux. Thus, exact solutions with energy flux can 
be now considered using the integrated framework. The fluid flow was calculated in terms of the metric functions 
for the double null, null, and diagonal coordinate systems. Then, the condition for perfect fluid was derived and 
other consistency relations were provided for each canonical type of coordinates. This work provides the needed part 
for an integrated framework to study and verify exact solutions with spacetimes of type B\ warped product. The 
approach developed is algorithmic and suited to the study of exact solutions using computer algebra systems such 
as the Interactive General Relativity Geometric Database [27J]. The framework developed is not limited to comoving 
reference frames. 
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VIII. APPENDIX: EINSTEIN AND WEYL TENSOR STRUCTURE FOR WARPED PRODUCT 

SPACETIMES OF CLASS Bi 

A. Double-Null Coordinates 

The nonvanishing components of the Einstein tensor are: 

_ 2(M,n B-B,, M,Q 
Gl1 - BM ' (92) 

_ 2MW 4 M : i2 —W 2 ,4 B + W,44 BW — W 2 ,s B + W,33 BW + 2M,i W 4 M,2 
12 ~ 2M 2 W i ' (93) 

= _2 { M,,B-B,M,) 

22 BM > V ) 

and 

„ (—B,i B,2 M + B,i2 BM + 2M, 12 B 2 )MW 2 

& 33 = G,44 — . (95) 

The nonvanishing components of the Weyl tensor components are all related and given by 

W 2 M 2 -B 

Cl 3 2 3 = ^1 4 2 4 = — Cl 2 1 2 = ^ Ci 4 3 4 (96) 

with 

Cl 3 2 3 = 

M 2 W A B,i2 B - M 2 W 4 B,i B,2 -2MW 4 M,i2 B 2 - B 3 W 2 ,* +B 3 W,*i W - B 3 W 2 ,a +B 3 W, 3 s W + 2B 2 M,i W i M,2 

6B 2 W 2 

(97) 

The condition for conformal flatness of the spacetimes is thus given by the vanishing of the numerator in (|9T[) . 
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B. Null Coordinates 

The nonvanishing components of the Einstein tensor are: 

G n = ; ~ 1 , J 2MW A M,n B 3 - 2MW i M,i B 2 B,i +MW 4 M,i B 2 A,2 -MW 4 M,2 B 2 A,i (98) 
M z W^B i 

+2MW i M,2 BAB,i +MW i M,2 BAA,2 -4AMW 4 M,i2 B 2 + 2MW 4 BA 2 M,22 -2MW 4 B, 2 A 2 M,2 
+A(W,* ) 2 B 3 - AW,** B 3 W + A(W,a ) 2 B 3 - AW,™ B 3 W - 2AM,i W 4 M,2 B 2 + M,l W 4 BA 2 ), 

G v , = (2MW i M,i2 B 2 - MW 4 M,2 A,2 B - 2MW 4 BAM,22 +2MW 4 B,2 AM,2 (99) 

-(W,4 ) 2 B 3 + W,^ B 3 W - (W,3 ) 2 B 3 + W, 33 B 3 W + 2M,i W 4 M,2 B 2 - M,1 W 4 BA), 

-2(M,22 B - B,2 M, 2 ) , 
G 22 = — y — - '—t, 100 

22 BM , V ) 



and 



MW 2 

G 33 = -^g-(4M,i2 B - 2A,2 M,2 B - 2B AM, 22 +2B,2 AM, 2 +2B MB, 12 -BMA,™ -25,2 MB,. +B,2 MA,2 



(101) 

The nonvanishing components of the Weyl tensor components are all related and given by 

W 2 M 2 B B —B 

Cl 3 2 3 = 2B Cl212 = -J Cl313 = ~7 Clili =Cn2i = 2 ]\ / pW 2 C3 4 3 4 ^ 102 ' 

with 



C1323 = ~^T7^oT77o (~2W, 2 B 3 + 2W,a B 3 W -2W, 2 B 3 + 2W, 33 B 3 W 
12B z \v 1 

-2(M,2 ) 2 W 4 BA - AMW 4 M,i2 B 2 + 2MW i M,2 A, 2 B + 4M,i W 4 M,2 B 2 + 2MFT 4 BAM,22 

-2MVK 4 B,2 AM, 2 2M 2 W i BB,i2 -M 2 W 4 BA, 22 -2M 2 W 4 B,2 B,i +M 2 W 4 B,2 A,2 ) (103) 

The condition for conformal flatness of the spacetimes is thus given by the vanishing of the numerator in (|103[) . 

C. Diagonal Coordinates 

The nonvanishing components of the Einstein tensor are: 

G n = tt^tttttt* (-MW i M,2 C,2 A + 2MW i M,22 AC + MW 4 M,i C,i C - (W,* ) 2 AC 2 (104) 
M z W^C Z 

+W,44 AC 2 W - (W, 3 ) 2 AC 2 + W, 33 AC 2 W + (M,i ) 2 W 4 C 2 + (Af, a ) 2 W i AC), 

_ (2M,i2 AC - A,2 M,i C - C,i M,2 A) 
Gl2 - MAC ' (105) 

G 22 = M 2^ 4A2 (~MW 4 M,i A,i C + 2MW 4 M, AC + MW 4 M,2 A,2 A - (W,4 ) 2 A 2 C (106) 



and 



+PT,44 A 2 CW - (W,s ) 2 A 2 C + W, 3 3 A 2 CW + (Af,i ) 2 W 4 AC + (Af, a ) 2 W 4 A 2 ), 



— MW 2 

G 33 = __(2A,i M,i C 2 - 4M )V C 2 A - 2A,2 M,a CM + 2C,2 M, 2 A 2 - 4M, 22 A 2 C (107) 
2(7,i M,i AC + (C,i ) 2 AM + (A, 2 ) 2 CM + A, 2 C,2 AM - 2A,22 CAM - 2C,n CAM + A,i C,i CM). 
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The nonvanishing components of the Weyl tensor components are all related and given by 

-W 2 M 2 _A A -A 

U 13 13 U 12 12 U 1414 ^^2 3 2 3 p U 2424 2W^ \ / 

with 

Ci 3 1 3 = 2AM l 2W 2 (4(W,a ) 2 A 2 C 2 - 4VK,44 A 2 C 2 VK - 414^,33 A 2 C 2 iy - 4(M,i ) 2 W 4 AC 2 - 4(M, 2 ) 2 irU 2 C 

+M 2 iy 4 (Ca ) 2 A + Af 2 iy 4 (A,2 ) 2 C + 4(Vy,4 ) 2 A 2 C 2 + 2AiW 4 M, 2 v4, 2 v4C* + 2AfW 4 Af,i C,i v4C* - 2MW 4 M,i A,i C 2 
+4MW 4 M,n AC 2 - 2MW 4 M,2 C, 2 A 2 + 4MW 4 M, 22 A 2 C + M 2 W 4 A,2 C, 3 A - 2M 2 W 4 A,22 CA - 2M 2 W 4 C,n CA + 
M 2 W 4 A,i C,i C) (109) 

The condition for conformal flatness of the spacetimes is thus given by the vanishing of the numerator in (|109| . 
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